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^ . ABSTRACT 

^ ' We use numerical simulations to investigate the origin and structure of the luminous halos 

that surround isolated galaxies. These stellar structures extend out to several hundred kpc 
' away from a galaxy, and consist of stars shed by merging subunits during the many accretion 

CS| \ events that characterize the hierarchical assembly of galaxies. Such origin suggests that outer 

luminous halos are ubiquitous and that they should appear as an excess of light over extrapo- 
lations of the galaxy's inner profile beyond its traditional luminous radius. The mass profile of 
\ the accreted stellar component is well approximated by a model where the logarithmic slope 

i steepens monotonically with radius; from p cx r^^ at the luminous edge of the galaxy to r^^ 

or steeper near the virial radius of the system. Such spatial distribution is consistent with that 
of Galactic and M3 1 globular clusters, suggesting that many of the globulars were brought 
■ in by accretion events, in a manner akin to the classic Searle-Zinn scenario. Luminous ha- 

ly-j I Ids are similar in shape to their dark matter counterparts, which are only mildly triaxial and 

f"^ . much rounder than dark halos formed in simulations that do not include a dissipative lumi- 

nous component. The outer stellar spheroid is supported by a velocity dispersion tensor with 
a substantial and radially increasing radial anisotropy; from cr^./cr^ 2 at the edge of the 
central galaxy to 5 at the virial radius. These properties distinguish the stellar halo from 
O ' the dark matter component, which is more isotropic in velocity space, as well as from some 

tracers of the outer spheroid such as satellite galaxies. Most stars in the outer halo formed 
C/3 ■ in progenitors that have since merged with the central galaxy; very few stars in the halo are 

^ ' contributed by satellites that survive as self-bound entities at the present. These features are 

in reasonable agreement with recent observations of the outer halo of the Milky Way, of M3 1 , 
, . and of other isolated spirals, and suggest that all of these systems underwent an early period 

' of active merging, as envisioned in hierarchical models of galaxy formation. 

b : 

5t 1 Key words: Galaxy: disk - Galaxy: formation - Galaxy: kinematics and dynamics - Galaxy: 

structure 



1 INTRODUCTION 

Galaxies have no edge. With rare exceptions, the stellar spatial dis- 
tribution in normal galaxies shows little sign of a sharp outer cutoff, 
and is reasonably well approximated by density laws that extend 
smoothly to arbitrarily large radius. Extrapolations of the inner lu- 
minosity profile, however, suggest that little light comes from re- 
gions of surface brightness much fainter than those traditionally 
used to define the luminous radii of galaxies (~ 25 mag/asec^). 
Perhaps because of this reason, together with the obvious obser- 
vational difficulties inherent to studying regions of low-surface 
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brightness, the outer luminous halos of external galaxies have in 
the past been regarded as a topic of little more than academic inter- 
est. 

This impression, however, is rapidly changing, as new datasets 
start to unveil some unexpected properties of the stellar component 
that populates the outer confines of galaxies. These developments 
have been made possible by the development of panoramic digi- 
tal cameras able to map the light distribution of external galaxies 
down to unprecedented surface brightness levels, complemented 
by efficient observational techniques designed to identify and mea- 
sure radial velocities of outer halo tracers in external galaxies, such 
as planetary nebulae (PNe, Romanowsky et al. 2003). Finally, the 
advent of efficient spectrographs in lOm-class telescopes have en- 
abled the measurement of radial velocities of large samples of giant 
stars throughout the Local Group (Ibata et al 2004), and dedicated 
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spectroscopic campaigns targeted on giant stars have dramatically 
increased the sample of tracers in the outer halo of our own Milky 
Way (Morrison et al 2003, Battaglia et al 2005, Clewley et al 2005). 

Likewise, the use of stacking techniques in wide-field surveys 
such as the Sloan Digital Sky Survey (SDSS) have led to clear de- 
tections of a distinct low surface brightness component around both 
isolated spiral galaxies (Zibetti, White & Brinkmann 2004) as well 
as central cluster galaxies (Zibetti et al 2005). The surface bright- 
ness profile of these outer luminous halos deviate significantly from 
straightforward extrapolations of the laws that describe the main 
body of the galaxy, and suggest that estimates of the amount of 
light in the intergalactic or intracluster medium might need to be 
revised upward. In the case of clusters, for example, Zibetti et al 
(2005) conclude that of order 11% of the light in a typical Abell 
cluster might be in the form of an intracluster luminous component 
with structural properties distinct from those of the galaxy cluster 
population and of the central cluster galaxy. 

These studies have brought about several unexpected results 
that challenge some of the accepted premises of the currently ac- 
cepted galaxy formation paradigm. For example, the kinematics of 
PNe in the vicinity of several normal ellipticals has been found to 
be consistent with simple models where such galaxies have no sur- 
rounding dark matter halo (Romanowsky et al. 2003, but see Dekel 
et al 2005). This result, if confirmed by further studies, would be 
very difficult to accommodate in the ACDM scenario, where all 
galaxies are envisioned to form as a result of the dissipative col- 
lapse of baryons within massive dark matter halos. 

In the case of the Local Group, wide field photometric surveys 
(and their follow-up spectroscopic campaigns) have modified radi- 
cally the traditional view of the outer stellar spheroid of the Milky 
Way and M31. Giant structures interpreted as tidal streams of dis- 
rupted or disrupting dwarfs have been identified in both galaxies 
(Ibata et al 1994, Helmi et al 1999, Majewski et al 2003, Ibata 
et al 2004), and have solidified the notion that accretion events 
may play a substantial role in the shaping of individual galaxies. 
These streams are expected to be short-lived as readily identifiable 
structures on the sky (Johnston et al 1995, Helmi & White 1999), 
and, although impressive, they are thought to contribute a relatively 
small fraction of stars in the halo of these systems. In broad terms, 
the main characteristics of the outer stellar halo are probably ad- 
equately captured by a simple model of a reasonably well-mixed 
spheroid of stars. 

Extensive spectroscopic surveys of giant stars in the outskirts 
of Local Group galaxies have unraveled some intriguing dynamical 
properties for the outer halo component. For example, just like PNe 
around ellipticals, the velocity dispersion of these halo tracers drops 
significantly in the outer regions of the Milky Way (Battaglia et al 
2005, Clewley et al 2005). This is a somewhat unexpected result if 
these galaxies are embedded in extended dark matter halos, and has 
prompted renewed interest in the extent and spatial distribution of 
dark matter in the outskirts of galaxies. 

From a theoretical point of view, the interest elicited by outer 
luminous halos stems from the fact that, in the currently accepted 
paradigm, stars are envisioned to form only in the collapsed, high- 
density regions near the center of dark halos — which we identify 
with the main body of individual galaxies. This prejudice is sup- 
ported by observations that indicate the need for a threshold gas 
density below which stars do not form (Kennicutt 1989, Martin & 
Kennicutt 2001, Schaye 2004), and imply that outer halo stars did 
not form in situ but have been shed from protogalaxies during the 
merger events that characterize the assembly of galaxies in a hier- 
archically clustering universe. 



In simple words, stars found as far away as, say, 100 kpc from 
a galaxy's center originate in satellites whose orbital apocenter was 
about that large before they spiraled in to merge with the main 
galaxy. We therefore expect a clear connection between the orbital 
properties of stars in the outer halo and those of the progenitors 
that merged to form the present-day galaxy. Unraveling the struc- 
ture of the outer luminous halo of galaxies may thus be seen as an 
important step toward unraveling the merging history of individual 
galaxies. 

In this paper, we analyze the origin and structure of the lumi- 
nous halos of galaxies simulated in the ACDM scenario. We de- 
scribe briefly in Section|2|the numerical simulations. The main re- 
sults of our analysis are presented in Section|3|and discussed in §|4] 
We conclude with a brief summary in Section|5| 



2 THE NUMERICAL SIMULATIONS 

We analyze a suite of eight numerical simulations of the formation 
of galaxies in the ACDM scenario. These simulations have been an- 
alyzed in earlier papers, which may be consulted for details on the 
code used as well as on the numerical setup (Steinmetz & Navarro 
2002, Abadi et al 2003a,b, Meza et al 2003, 2005). In brief, each 
simulation follows the evolution in a ACDM universe of a small 
region surrounding a target galaxy, excised from a large periodic 
box and re-simulated at higher resolution preserving the tidal fields 
from the whole box. The simulation includes the gravitational ef- 
fects of dark matter, gas and stars, and follows the hydrodynamical 
evolution of the gaseous component using the Smooth Particle Hy- 
drodynamics (SPH) technique (Steinmetz 1996). We adopt the fol- 
lowing cosmological parameters for the ACDM scenario: Ho ~ 65 
km/s/Mpc, erg = 0.9, = 0.7, ^cdm = 0.255, fibar = 0.045, 
with no tilt in the primordial power spectrum. All simulations start 
at redshift Zinit — 50, have force resolution of order 1 kpc, and 
mass resolution so that each galaxy is represented, at z — 0, with 
~ 125, 000 star particles. 

Dense, cold gas in collapsed regions is allowed to turn into 
stars at rates consistent with the empirical Schmidt-like law of Ken- 
nicutt (1998). The energetic feedback of evolving stars is included 
mainly as a heating term on the surrounding gas, but the short 
cooling times in these regions reduce significantly the effective- 
ness of feedback in curtailing star formation. The transformation 
of gas into stars thus tracks closely the rate at which gas cools and 
condenses at the center of dark matter halos. As discussed in the 
references listed in the above paragraph, this results in an early on- 
set of star-forming activity in the many progenitors of the present- 
day galaxy present at high redshift. This leads to the formation of 
a prominent spheroidal component as these progenitors coalesce 
and merge to form the final galaxy. Gas accreted after the merging 
activity abates leads to the formation of a centrifugally-supported 
gaseous and stellar disk component clearly present in most of our 
simulations, although they make up typically only about ~ 30% of 
the final stellar mass of each galaxy. 

It is likely that improvements to our feedback algorithms may 
lead to revisions in the efficiency and timing of star formation 
in these galaxies. As we discuss below, outer halo stars originate 
mainly in the merger of the early progenitors, so reducing the ef- 
ficiency of early star formation should have some impact on the 
actual prominence of the outer stellar halo component. We expect, 
however, that reasonable modifications of the star formation algo- 
rithm will affect mainly the number, ages, and chemical compo- 
sition of outer halo stars, rather than its dynamical and structural 
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Figure 1. Spatial distribution of the dark matter and stellar component of one simulated galaxy at z=0. Top panels show the dark matter (left) and the stars 
(right) within the virial radius of the system. Dark matter particles are colored by their local density, whereas stars are colored by age (blue: 0.0 < age/Gyr < 
2.5, cyan: 2.5 < age/Gyr < 5.0, green: 5.0 < age/Gyr < 7.5, yellow: 7.5 < age/Gyr < 10.0, red 10.0 < age/Gyr < 15.0). The outer green circle shows the 
virial radius of the system, the inner circle shows the radius adopted as the "luminous radius" of the galaxy. Bottom panels split the stellar component into two 
different groups: stars formed in the inain progenitor of the galaxy (i.e., "in situ" stars, left panel) and those contributed by accretion events ("accreted" stars, 
right panel). The latter group excludes stars in satellites that remain self-bound in the halo of the galaxy. 

properties, which we expect to depend mainly on the orbital prop- lar speeds comparable to the Milky Way. We summarize the main 
erties of the merging progenitors and on the structure of the dark structural parameters of the dark matter and stellar components of 
matter halo host. These properties are less sensitive to the complex the simulated galaxies in TableQ 
astrophysics of star formation and feedback, and we therefore fo- 
cus our analysis on the structure and dynamics of the outer stellar 
halo in our eight simulations. These target a small range in mass 
chosen so that sA z — the galaxies have luminosities and circu- 
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Figure 2. The radial dependence of the fraction of stars formed in the most 
massive progenitor of the final galaxy (i.e., stars formed in situ) is indi- 
cated by dotted lines connecting symbols. The thick solid lines indicate the 
cumulative fraction of stars within the virial radius. In-situ stars dominate 
throughout the main body of the galaxy; within the "luminous radius" of the 
galaxy, defined to be rju^ = 20 kpc (shown with arrows), they contribute 
~ 60% of the stellar mass. Outside ri^u the stellar component consists 
almost exclusively of accreted stars; indeed, fewer than ~ 5% of stars in 
the outer (r > r'lum) halo formed in situ. 



3 RESULTS 

Figure Q shows, at 2 = 0, one of our simulated galaxies (KIA3, 
see Table projected onto a box of 540 kpc on a side. The top 
panels show the dark matter particles (left) and stars (right) within 
the virial radius (r^ii ~ 270 kpc, shown by the outer green circle), 
defined to encompass a region of mean density 100 times the criti- 
cal density for closure. Dark matter particles are coloured by their 
local density, while stars are coloured by their age, as described in 
the label. 

The bottom panels in FigureQseparate the stars in two com- 
ponents: "in situ" stars that formed in the most massive progenitor 
(left) and "accreted stars" that formed in progenitors that merged 
with the main galaxy (right). Stars labelled as "accreted" exclude 
those associated with self-bound satellite systems that survive un- 
til the present. These can be seen clearly in the top-right panel of 
Figure Q but are largely absent in the bottom-right panel, except 
for a sprinkling of stars associated with a young stream recently 
torn from a disrupting satellite, visible at a; ~ —20, j/ ~ 10 kpc. 
Roughly ~ 48% of stars (by mass) formed "in situ" in this galaxy, 
compared with ~ 44% which make up the accreted component. 
Satellites contribute a rather small fraction (~ 8%) of all stars 
within Tvir- These numbers are typical of our simulations, as may 
be seen from the numbers listed in TableQ 

FigureQillustrates a number of properties of the stellar com- 
ponent common to all of our simulations. In particular, it is impor- 
tant to note that (i) stars spread as far out as the virial radius of 
the system (outer circle in the panels of Figure 0, although they 
are more highly concentrated than the dark matter halo; (ii) in situ 



Figure 3. Density profile of stars (circles) and dark matter (squares) in one 
of our simulations (KIA3), shown in a logarithmic scale. Open circles cor- 
respond to the stars formed in situ, whereas filled circles correspond to ac- 
creted stars. Note that the shape of the density profile of the accreted stellar 
component is similar to that of the dark matter, and can be adequately fit 
with a radial law where the logarithmic slope is a power law of radius. See 
text for details. 




Log(R/kpc) 



Figure 4. i?-band surface brightness profile of stars, split between the "in- 
situ" and "accreted" components (run KIA3, see Table IT]. The accreted 
stars dominate outside the luminous radius and appear as an excess of hght 
over a bulge-l-disk fit to the inner surface brightness profile, which is shown 
as a solid line. The dashed line shows a Sersic-law fit to the outer profile 
of accreted stars. Note that a single Sersic-law (with parameters given in 
Table|5| reproduces very well the radial distribution of all accreted stars. 
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stars (bottom left panel), are responsible for most of the young stars 
in the main body of the galaxy (note the prominence of the young 
disk in the bottom-left panel), and are practically absent from the 
outer halo; whereas (iii) accreted stars (bottom right panel) make 
up preferentially the spheroidal component and dominate the stel- 
lar budget in the outer regions of the galaxy. 

The relative importance of the "accreted" vs "in situ" com- 
ponents at various radii is shown in Figure |2| Beyond a radius 
''lum ~ 20 kpc (marked by arrows in Figure|2j accreted stars dom- 
inate in all our simulations. We shall define rium as the "luminous 
radius" of the galaxy and refer to stars beyond rium as the "outer 
luminous halo" or "outer galaxy", for short. Our analysis focuses 
on the structure and dynamics of the outer halo and, consequently, 
mainly on the properties of the accreted component. 

It is also clear from this figure that the outer luminous halo 
contains a relatively small fraction of the stellar mass of the galaxy; 
fewer than 15% of all stars are found between rium and the virial 
radius. We note as well that accreted stars make up a non-negligible 
fraction of stars in the inner galaxy (i.e., inside rium); we shall 
therefore distinguish between "accreted" and "outer halo" stars in 
what follows. 

Figure|3|shows the density profile of the stellar and dark mat- 
ter components of the simulation shown in Figure Q The stellar 
contribution is split between in-situ and accreted stars. Note that 
stars are significantly more concentrated than the dark matter; in- 
deed, they contribute a significant fraction of the total mass within 
the luminous radius (44% on average for all runs in our series), but 
they make up a negligible amount of the total mass in the outer re- 
gions. Note as well the sharp truncation of in-situ stars at r > rium; 
essentially no stars in the outer halo were formed in the main pro- 
genitor. This is a natural consequence of the fact that in-situ star for- 
mation proceeds efficiently only in regions of high-density. The few 
in-situ stars found beyond rium at 2: = have been torn from the 
main progenitor during mergers, although Figure|2|(and the data in 
Table0 shows that this is rather inefficient, and that the "accreted" 
component dominates the stellar component of the outer halo. 

The radius beyond which accreted stars begin to dominate is 
signaled by an abrupt change in the surface brightness profile of 
the galaxy, as shown in Figure|4| Inside the luminous radius, where 
fiR is brighter than about 25.5 mag/asec^, the surface brightness 
profile of this galaxy (projected in this case face-on and studied in 
detail by Abadi et al 2003a,b) is well approximated by a de Vau- 
couleurs bulge plus an exponential disk (shown as a solid line). 
This bulge-l-disk model is, however, unable to fit the structure of 
the outer luminous halo, which deviates from the inner profile and 
appears as a outer luminous "excess" beyond ~ 30 kpc, where the 
7?-band surface brightness drops to values fainter than about 26.5 
mag/asec^ . 

The outer halo surface brightness profile steepens gradually 
outwards and is well approximated by a Sersic-law, 

E(i?) = E,ffe-^"'(«/«"'-f''"'-^l (I) 

where Rett and Ecfi are, respectively, the radius containing half 
the light and the surface brightness level at that radius. The surface 
brightness profile of the outer regions of our simulated galaxies 
is shown in Figure [S^, together with Sersic-law fits to the outer 
portion of the profiles. 

The Sersic fit parameters are listed in Table|2| On average we 
find (n) = 6.3; (i?off) = 7.7 kpc. The outer luminous halo's sur- 
face brightness profile steepens from E(i?) oc R~^'^ at about the 
luminous radius to R~^'^ at r ~ 100 kpc and to R~^'"' or steeper 
around the virial radius. The gradually steepening profile of the 



outer halo characterized by the Sersic law is reminiscent of the dis- 
tribution of dark matter, whose density profile also steepens mono- 
tonically outward (Navarro et al 2004). The dark matter profile is, 
however, much less centrally concentrated and its slope is shallower 
than the stellar halo at all radii. As discussed by Merritt et al (2005), 
the dark matter density profile may also be approximated by a Ser- 
sic law, but with a characteristic value of 2 < n < 4 and much 
larger effective radii. 

Figure|5l' shows the surface brightness profiles of the accreted 
stellar component only for all eight simulations in our series. The 
profiles have been scaled so that they would all line up along a line 
of slope -1 if they followed accurately a Sersic law. All profiles are 
indeed well approximated by eq. 1 , with the possible exception of 
one case where a "bump" of stars associated with an ongoing satel- 
lite disruption event is seen just outside R^s (open starred symbols 
in Figure|5j>). This raises the interesting prospect of estimating the 
total fraction of stars accreted throughout the history of a galaxy 
simply by fitting its outer profile. As the last two columns of Ta- 
ble |2| demonstrate, the total light contributed by accreted stars is 
well approximated (to within a factor of ~ 2) by the total light of a 
Sersic-law fit to the outer luminous halo. 

The similarity between profiles suggest that the violent relax- 
ation associated with mergers endows the accreted stellar compo- 
nent with a simple, roughly self-similar structure that is well ap- 
proximated by a Sersic law. This process has been studied exten- 
sively in the literature, and previous studies have consistently found 
that the structure of N-body merger remnants is quite reminiscent 
to that of bright ellipticals (see, e.g., Barnes & Hemquist 1992 and 
references therein), which are indeed well approximated by a Sersic 
law with n>4 (see, e.g., Graham & Guzman 2003 and references 
therein). 

One may think of the accreted stellar component as formed by 
the overlap of the many "tidal tails" stripped from the merging pro- 
genitors; in this interpretation the outer luminous halo would just 
be the superposition of the outer tails stripped from each progeni- 
tor during the merger process. This characterization of the accreted 
stellar component helps to explain the remarkable kinematics of the 
outer halo shown in Figure|6| This figure shows the velocity disper- 
sion profile of stars and dark matter, averaged over all simulations 
after scaling the positions and velocities of all particles to the virial 
radius and virial velocity of each system. The right panels in Fig- 
ure|6|show that the velocity dispersion of the dark matter declines 
gradually from the center outwards, and is characterized by a mild 
radial anisotropy, going from almost isotropic near the center to 
P — {1 — (T^/(T^) ~ 0.3 in the outer regions. (We use at to denote 
the tangential velocity dispersion; = (erg + cr^)/2.) 

The outer luminous halo, on the other hand, shows a much 
more pronounced radial bias: from P ~ 0.4 just outside the lumi- 
nous radius of the galaxy (rium is on average of order 0.06 rvir) to 
about /3 ^ 0.8 near the virial radius. This substantial-and mono- 
tonically increasing-radial anisotropy may also be understood as 
a direct consequence of the tidal disruption process of formation 
described above. Outer halo stars are typically stripped from merg- 
ing progenitors during pericentric passages. However, because stars 
form in high-density regions the stripping of stars operates effi- 
ciently only when pericenters are small (and tidal effects are great- 
est). This occurs typically after dynamical friction has eroded the 
orbit sufficiently to bring the pericenter of the accreting satel- 
lite close to the central galaxy. As a result, outer halo stars are 
"launched" into highly energetic orbits with, characteristically, the 
(small) pericentric radii that accompany the most disruptive tides. 
Stars able to reach farther are thus typically on more eccentric or- 
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Figure S. (a-left) R-band surface brightness profile of the outer galaxy (10 < -R/kpc< 300) in all of our simulations. Solid lines correspond to the best 
Sersic-law fits to the outer regions of the profile (r > rium) whose parameters are listed in Table|2] (b-right) Same as in the left, but profiles correspond only 
to accreted stars and have been scaled as stated in the axis labels so that Sersic-law profiles would follow the straight solid line of slope —1. Interestingly, 
a simple Sersic-law reproduces quite well the radial distribution of accreted stars in all simulations. There is one exception (which shows as a bump in the 
right-hand panel), associated with the transient effect of an ongoing satellite disruption event. 



bits, leading to the increasing radial anisotropy in the orbits of the 
outer halo stars seen in Figure|6| 

The formation process of the accreted stellar component is 
thus qualitatively similar to that of the dark matter halo, where 
mergers also play a substantial role. It is thus perhaps not sur- 
prising that they settle onto similar, mildly triaxial structures, as 
shown in FigureQ For the accreted stars, the average intermediate- 
to-major axis ratio is (b/a) = 0.91, whereas the average minor- 
to-intermediate axis ratio is (c/&) = 0.92 ({c/a) = 0.84). This 
is not too different from the dark matter halos, which have on av- 
erage {b/a) = 0.94 and (c/b) = 0.90 {{c/a) = 0.84). We note 
that the dark matter halos in these simulations are much rounder 
than typically found in N-body simulations of CDM halo forma- 
tion: (b/a) ~ 0.75 and (c/a) ~ 0.6 (see, for example, Bailin & 
Steinmetz 2004). This is the result of the response of the halo to the 
dissipative assembly of the baryonic component of the galaxy at its 
center, which steepens the potential well and reduces the triaxiality 
of the halo; see, e.g., Kazantzidis et al (2004), Bailin et al (2005), 
and references therein. 

The in-situ stellar component, on the other hand, often sports a 
well defined rotationally-supported disk (see, e.g., lower left panel 
of Figure0 and its shape is thus well described by an oblate struc- 
ture with (b/a) — 0.95 and (c/a) — 0.62. Although rotation plays 
little role supporting the accreted component (the mean rotation ve- 
locity, = jz/R, is typically less than one tenth of the circular 
velocity at each radius), it is interesting to note that the direction of 
its angular momentum is well aligned with that of the inner galaxy. 

This is shown in Figure|8| which show the distribution of the 
cosine of the angles between the angular momentum of the dark 
matter halo (Jdrk), the in-situ stars ( Jins), and the accreted stellar 
halo (Jacc)- Although the good agreement between Jacc and Jins 
could have perhaps been anticipated (after all, they are both part 



of the same stellar system) it is still interesting to note the strong 
correlation between the rotational properties of the iimer galaxy 
and the dark matter halo. This implies that, although mergers lead 
to substantial transfer of angular momentum from the stars to the 
dark halo (Navarro & Steinmetz 1997, Abadi et al 2003a), this does 
not alter radically its orientation, so that the spin of the inner galaxy 
generally ends up aligned with the rotation axis of the surrounding 
dark matter halo. 

It is also of interest to characterize the orientation of the rota- 
tion axis relative to the principal axes of the system. Figure|9|shows 
the alignment between the galaxy's angular momentum and the mi- 
nor axis of the dark matter and luminous halos. No obvious align- 
ment is seen in these panels, suggesting that the rotational proper- 
ties of the galaxy are, at best, weakly correlated with the shape of 
the surrounding dark matter (or luminous) halo. We shall return to 
possible interpretations of this result in the following section. 

We end the characterization of the outer luminous halos of 
simulated galaxies by comparing the age distribution of stars in the 
inner and outer galaxy, as well as that of satellites enclosed within 
the virial radius (Figure fTol . The ages of stars in the outer galaxy 
(r > rium) differ significantly from those in higher density regions 
(like the inner galaxy or the surviving satellites), where star forma- 
tion may proceed. The outer halo is populated mainly by older stars, 
reflecting the fact that the mergers responsible for its formation are 
more common at earlier times. Interestingly, the distribution of ages 
of stars in the outer halo is also fairly distinct from that of stars in 
satellites orbiting within the virial radius. This shows that relatively 
few stars in the outer halo originate in the "harassment" of satellites 
that have survived as self-bound entities until the present (Moore et 
al 1996). Most stars in the outer halo come from merger events 
whose progenitors have long been fully disrupted, suggesting that 
the properties of the satellite population may be quite distinct from 
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Figure 6. Average velocity dispersion profile of stars and dark matter in our 
simulations. Averages are computed by scaling all positions and velocities 
to the virial radius and virial circular velocity of each system, respectively. 
The rotation axis of the inner galaxy is chosen as the polar axis of the scal- 
ing procedure. Particles are then grouped in spherical bins and the spherical 
components of the velocity dispersion are computed. These profiles are fi- 
nally averaged over all systems to produce the radial and tangential velocity 
dispersion profiles plotted here. Right panels con'espond to the dark mat- 
ter component and left panels to the accreted stars. Recall that, on average, 
the luminous radius is ~ 0.06 rvir in these units. Note that the dark mat- 
ter velocity dispersion exhibits a mild radial bias; the radial bias is much 
more pronounced, and monotonically increasing, in the case of the accreted 
stellar component. 

that of the smooth outer halo. We explore the consequences of these 
results for the interpretation of observational data in the following 
section. 



4 DISCUSSION 

The results presented in the previous section have a number of im- 
plications regarding the interpretation of observations of stars in re- 
gions far removed from the normal boundaries of typical galaxies. 
Our simulations suggest that these "intergalactic stars" (or, more 
properly, "extra-galactic stars") are a relict of the merging history 
of each individual galaxy. As such, all galaxies with a past of active 
merging activity (the majority in a hierarchically clustering sce- 
nario such as ACDM) should have an outer stellar halo component 
with properties similar to those described in the previous section. 

The outer luminous halo should appear as an excess of light 
over the extrapolated profile of the inner galaxy. Its structure is 
characterized by a gradually steepening density profile and by a 
clearly defined radial anisotropy in the velocity distribution. Outer 
halos are rather difficult to detect, since they are-by definition- 
confined to regions of very low surface brightness difficult to probe 
observationally. However, the presence of outer luminous compo- 
nents detached from the properties of the inner galaxy have been 
recently reported by a number of authors. Following the pioneering 
work of Sackett et al (1994), Morrison et al (1994), and Zheng et 



Figure 7. Inertia axis ratios of the dark matter halo, as well as of the accreted 
and in-situ stellar components. Axis ratios are computed by diagonalizing 
the inertia tensor lij = J] rUiXiXj for all pailicles of each component 
within the virial radius of the system. Note that the shape of the dark mat- 
ter and accreted components are similar and only mildly triaxial, but that 
the in-situ stars are predominantly oblate. This reflects the fact that in-situ 
stars have a higher proportion of newly formed stars, which tend to arrange 
themselves in a disk-like structure (see FigurefT]. 

al (1999) among others, Zibetti, White & Brinkmann (2004), have 
recently reported the statistical detection of a p oc r^^ halo of stars 
in the outskirts of a sample of edge-on disk-dominated galaxies se- 
lected from the Sloan Digital Sky Survey (SDSS). 

Encouragingly, this outer halo component is detected as ex- 
cess light over an extrapolation of the inner disk (see Figure|4}; the 
halo dominates the minor axis light profile beyond about ~ 6 ex- 
ponential scalelengths, in regions where the r-band surface bright- 
ness drops to values fainter than fir ~ 27.5. A qualitatively simi- 
lar result has been derived for M3 1 through detailed stellar counts 
by Guhathakurta et al (2005) and by Irwin et al (2005), who find 
that the surface brightness profile along the minor axis "flattens" 
relative to an extrapolation of the de Vaucouleurs fit to the inner 
spheroid at radii beyond ~ 20 kpc. 

Our simulations provide a compelling interpretation where 
outer stellar halos are made up predominantly of stars accreted dur- 
ing previous merger events and, in particular, by the subset which 
were propelled into highly energetic (and highly eccentric) orbits 
during the disruption process that accompanies the mergers. Merg- 
ers and phase mixing lead then to the formation of a tenuous, dis- 
tinct stellar component that fills the halo of the galaxy out to the 
virial radius. 

Although we expect them to be ubiquitous, the prominence of 
these outer luminous halos is difficult to assess from a theoretical 
point of view, as it depends critically on the number and timing of 
merger events, as well as on the fraction of stars present at the time 
the mergers occur, all of which remain highly uncertain. A number 
of general inferences are, however, still possible. Late major merg- 
ers should propel a proportionally larger fraction of stars into the 
outer halo, and therefore the relative prominence of this component 
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mentum vector of the in-situ stellar component, Jins, the accreted stars, 
Jacc: and the dark matter, Jjik- Each of the three different permutations 
are shown by different shaded histograms, as indicated in the labels of each 
panel. 



should be higher in spheroid-dominated galaxies, which are widely 
regarded as merger remnants. This process should be especially im- 
portant in systems that have undergone repeated mergers, leading 
to the expectation that the relative importance of the "intracluster" 
stellar component should be highest in clusters with prominent cD 
galaxies. 

In disk-dominated galaxies, where in-situ star formation is 
typically still ongoing and mergers are less important, the outer 
halo is likely to dominate the profile only at large radii, where star 
formation thresholds prevent the efficient formation of in-situ stars 
(Kennicutt 1989, Ferguson et al 1998, Schaye 2004). The accreted 
stellar component — to which no such threshold applies — appears 
then as a distinct outer component, "flattening" the surface bright- 
ness profile in the outer regions, as reported by Zibetti et al (2004), 
Guhathakurta et al (2005) and Irwin et al (2005). 

Quantitatively, we expect the slope of the surface brightness 
profile of the outer luminous halo to steepen gradually outwards, 
as dictated by the Sersic-law fits summarized in Table|2| Given that 
the average effective radii of the accreted component is 7.7 kpc and 
(n) = 6.3, we expect its surface brightness profile to steepen from 
from E oc R^^'^ at the luminous radius (20 kpc) to E oc R~^'^ 
at r ~ 100 kpc. This is in reasonable agreement with the results 
of Zibetti et al (2004), Irwin et al (2005), and Guhathakurta et al 
(2005). The latter authors, in particular, report that the M31 outer 
halo may be approximated by E oc R~^"^ in the range 30 kpc 
< R < 100 kpc, slightly shallower than in our simulations, but not 
inconsistent given the substantial observational uncertainty. 

We note as well that the distribution of globular clusters in the 
Milky Way and M31 is consistent with that of the accreted stel- 
lar component in the simulations. With the caveat of small num- 
ber statistics (only ~ 150 globular clusters around the Galaxy and 
~ 360 around M31 are currently confirmed; see Harris 1996 and 
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Figure 9. (a)-top: Distribution of the cosine of the angle between the an- 
gular momentum of in-situ stars, Jins, accreted stars, Jacc, and the dark 
matter, Jdrk. ™d the minor axis of the dark matter . As in Figure |8] his- 
tograms of different shades refer to the various permutations, as indicated 
in the labels, (b)-bottom: Same as (a), but for the minor axis of the accreted 
stellar component. 



Galleti et al 2004), it is interesting that a Sersic law with n ~ 6 is 
consistent with the (projected) number density profile of Galactic 
globulars. This is shown in Figure fTTI where the profile of M31 and 
Galactic globulars is compared with the Sersic-law that describes 
the average accreted stellar component in our simulations. We em- 
phasize that the curves shown are not fits, and therefore the general 
agreement is quite suggestive. A simple interpretation, of course, 
is that the globular cluster population shares the same "external" 
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origin as that of accreted stars in our simulations, and that they are 
thus relicts of the accretion history of each galaxy. (Guhathakurta 
et al (2005) note as well that the smooth halo of M3 1 red giant stars 
they detect out to ~ 150 kpc mimics the distribution of its globular 
cluster population.) Taken together, the evidence linking globular 
clusters with accretion events, as envisioned in the classic model of 
Searle & Zinn (1978), seems quite strong. 

In terms of space density, the slope of the accreted stellar com- 
ponent is close to p oc r^"*'^ near the edge of the luminous galaxy. 
This is in reasonably good agreement with that inferred from trac- 
ers of the stellar halo of the Milky Way in the vicinity of the Sun 
(see, e.g., Zinn 1985, Saha 1985, Morrison et al 2000). Our sim- 
ulations imply that the halo profile should steepen gradually with 
radius and become significantly steeper than r""^ in the outer re- 
gions. It is interesting to explore the consequences of this gradual 
steepening for the interpretation of dynamical data in the outskirts 
of the Milky Way. 

Battaglia et al (2005) report that the velocity dispersion of 
halo tracers such as blue horizontal branch stars and red giants de- 
clines from ~ 120 km/s at the edge of the luminous disk to ~ 50 
km/s at ~ 120 kpc. These authors argue that the data favour a dark 
matter halo model sharply truncated in the outer regions over the 
extended mass distributions predicted for CDM halos, such as the 
Navarro, Frenk & White (1996, 1997, NFW) profile. Neither halo 
model, however, is completely satisfactory. Extreme radial velocity 
anisotropy ((3 ~ 1) is required to reconcile the data with an NFW 
halo model although the opposite (i.e., a tangentially anisotropic 
velocity distribution, /3 ~ —0.5) is required for the truncated halo 
model. The latter seems an unlikely result given that the outer re- 
gions of systems assembled hierarchically are typically dominated 
by radial motions (see, e.g., Hansen & Moore 2004 and references 
therein). 

We can use our results to explore the consistency of the veloc- 
ity dispersion data with CDM halos. We write the spherical Jeans' 
equation as. 



dlnr 



(2) 



where we have used a * subscript to denote quantities associated 
with the stellar halo. The terms in the left-hand side of eq|2|may be 
evaluated by combining the Battaglia et al data with the assumption 
that the stellar halo profile may be approximated by a Sersic-law 
with the (average) parameters given in Table |2| The result of this 
exercise is shown in Fieure fT2l 

Interestingly, this simple model predicts a circular velocity 
near the center consistent with that measured at the solar circle 
(~ 220 km s~^) without tuning (solid squares in the bottom panel 
of Fieure ll2> . Further out, the sustained drop in cr* leads to much 
lower circular velocities, approaching ~ 150 km s^^ at r ~ 100 
kpc. This drop in the circular velocity of the system illustrates why 
Battaglia et al's analysis favours a truncated halo. We argue, how- 
ever, that the data are consistent with an extended CDM halo, and 
that the velocity dispersion drop places interesting constraints on 
its total mass. 

This is shown in the bottom panel of Fieure [T2l bv the top thick 
solid line, which corresponds to an (adiabatically contracted) NFW 
halo of virial velocity Vvir ~ 110 km/s and concentration c — 
14 (the average for halos of this mass in the ACDM cosmogony). 
The halo contribution before and after contraction is shown by the 
dotted and dashed lines, respectively. 

The relatively good agreement with the Vc profile derived 
from the data implies that, rather than being truncated, the halo 
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Figure 10. Age distribution of stars in the inner galaxy (i.e., stars within 
riuni, solid blue line), in the outer galaxy (r > rim,j, dotted red line) and 
in satellites (dashed green line) enclosed inside the virial radius r^i^ and 
averaged over our set of 8 simulations. 



of the Milky Way may just be less massive than commonly as- 
sumed. The virial mass of the fit shown in Figure fT2l is only ~ 
6.4 X 10^^ Mq, and its virial velocity is only one half of the circu- 
lar speed at the solar circle. This result is similar to that of Klypin, 
Zhao & Somerville (2002), who analyzed constraints placed by dy- 
namical tracers within the luminous radius of the Milky Way and 
found that reconciling such data with CDM halos requires the virial 
velocity of the dark matter halo to be substantially lower than the 
rotation speed of the disk. (Their favourite Milky Way halo model 
has Mvir ~ 10^^ Mq, which corresponds to r^h- ~ 260 kpc, and 
Vvir ~ 130 km s~^.) One may rephrase this result by stating that 
the relatively small dark matter contained within the luminous ra- 
dius of the Galaxy is easiest to accommodate with low-mass halo 
models, which have the virtue of predicting directly a sharp drop 
in the velocity dispersion of outer halo tracers. A similar argument 
was developed by Navarro & Steinmetz (2000), although note the 
erratum presented in Eke, Navarro & Steinmetz (2001). 

It is unclear whether the mass of the halo of the Milky Way 
is unusually low compared with other spiral galaxies of similar lu- 
minosity, but galaxy-galaxy weak lensing data seem to favour halo 
masses of the order we find here for late-type L* galaxies (see, e.g., 
Guzik & Seljak 2002). Also, there is now convincing evidence that 
the velocity dispersion of satellite galaxies drops at large projected 
radii so the drop in Vc in the outer regions may actually be a com- 
mon feature in Milky Way-like galaxies (see Brainerd 2004 for a 
review). 

Turning our attention now to the velocity structure of the outer 
halo, we recall the pronounced radial anisotropy shown in Figure|5| 
and argue that such effect should not be ignored when interpreting 
the dynamics of halo tracers in the outskirts of galaxies. As shown 
by Dekel et al (2005), accounting for such anisotropics may be 
enough to reconcile the steeply declining velocity dispersion pro- 
files of planetary nebulae around bright ellipticals (Romanowsky 
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Figure 11. Projected number density profile of globular clusters around the 
Milky Way (filled circles) and M31 (filled squares), respectively. Data for 
the 147 Galactic globulars and for the 363 M31 globulars are from Hanis 
(1996) and Galleti et al (2004). We assume that 100 arcmin= 22 kpc at the 
distance of M3 1 . We have averaged the distribution of Galactic globulars 
over several random lines of sight to construct the projected number den- 
sity profile. Each bin contains 15 and 40 globulars for the Galaxy and M31, 
respectively, and has Poisson en'or bars. The two curves are not fits to the 
data; rather, they illustrate the distribution of accreted stars in our simula- 
tions. Parameters for the two Sersic-law curves are chosen to be the average 
of the data presented in Tablel2land are vertically scaled to match the data 
for each galaxy. 



et al. 2003) with the presence of massive dark halos around these 
objects. The radial anisotropy differentiates dynamically the stars 
in the outer halo from the dark matter (Figure |6j. It also distin- 
guishes the stellar halo from the population of surviving satellites, 
which are found to trace rather faithfully the main properties of the 
dark halo. Indeed, the average anisotropy of the satellite population 
in our simulations is found to be /3 ~ 0.45 (Sales et al, in prepara- 
tion), closer to that of the dark matter, and certainly less anisotropic 
than the accreted stellar component. 

This dynamical distinction is not surprising, as we find that 
the smooth stellar halo has little relation with the population of sur- 
viving satellites (see Figure [Tol . Indeed, few stars in the outer halo 
may be traced to satellites that survive as self-bound entities until 
the present. The properties of the stellar halo seems more closely 
linked to the progenitors of early mergers than to the "harassment" 
of the satellite population that survives until the present day. 

This result has interesting consequences when applied to the 
interpretation of the origin of the intracluster light (ICL) compo- 
nent. Zibetti et al (2005), for example, find that the properties of 
the ICL in their sample of clusters selected from the SDSS are quite 
different from those of the cluster galaxy population as a whole. In 
particular, the ICL is found to be more centrally concentrated than 
the cluster galaxies and more significantly aligned with the (bright- 
est cluster galaxy (BCG) than cluster galaxies. In other words, the 
ICL seems more intimately related to the central galaxy (BCG) than 
to the cluster galaxy population. All of these properties are consis- 



tent with a scenario where the origin of the ICL is traced to the 
mergers that led to the formation of the BCG, with a minor contri- 
bution from stars stripped from surviving galaxies. 

A final item for discussion is the weak alignment between dark 
halo shapes and angular momentum shown in Figure|9| At first, this 
seems in disagreement with previous studies, many of which report 
significant alignment between the rotation axis and the minor axis 
of CDM halos (see, e.g., Bailin & Steinmetz 2004 and references 
therein). The issue may be resolved by noting that the shapes of the 
dark halos change dramatically when a dissipative baryonic compo- 
nent is included; indeed, in most cases the dark halos are so nearly 
spherical that the precise directions of the principal axes are poorly 
determined. Furthermore, as discussed by Bailin et al (2005), it ap- 
pears as if the inner regions of the dark matter halo respond to the 
assembly of the galaxy by aligning its minor axis with that of the 
luminous component, but at the same time decouple from the ori- 
entation of the halo in the outer regions. 

Interestingly, the weak correlations between halo shape and 
spin that our simulations predict could be tested directly with up- 
coming weak lensing surveys. Hoekstra et al (2004) report a tan- 
talizing correlation between the shape of the luminous galaxy and 
that of their surrounding matter, but their sample contains a broad 
mixture of morphological types and a relatively strong contribution 
from spheroid-dominated galaxies, where the relationship between 
shape and spin is less clear than in spiral galaxies. Weak lensing 
analysis applied to spiral-dominated samples should provide fur- 
ther constraints on this topic. This is currently being addressed 
using large, morphologically-segregated samples from the SDSS 
(Mandelbaum et al, in preparation). 



5 SUMMARY 

We investigate the origin of stars that populate regions well beyond 
the traditional luminous boundaries of normal galaxies using nu- 
merical simulations of galaxy formation in the ACDM cosmogony. 
Our simulations shows that such components are ubiquitous and 
owe their presence to the many mergers that characterize the early 
formation history of galaxies assembled hierarchically. These faint 
stellar halos extend out to the virial radius of the system, and con- 
sist mainly of the overlap and mixing of the many "tidal tails" shed 
by merging progenitors during the assembly of the galaxy. Such 
origin leads to robust predictions for the dynamics and structure of 
the outer spheroid that may be contrasted with observation. 

• The density profile of the accreted stellar component — which 
dominates the outer luminous halo — is well approximated by a 
Sersic-like model where the logarithmic slope steepens monoton- 
ically with radius; from p cc r^^ near the edge of the galaxy's 
traditional luminous boundary to r~* or steeper near the virial ra- 
dius. The shape and concentration of the accreted component is in 
reasonable agreement with the Galactic and M31 globular cluster 
population, lending support to the classic Searle-Zinn scenario for 
the formation of the globular cluster population. 

• The accreted stellar component is reasonably well approx- 
imated with a mildly triaxial spheroid with average axis ratios 
{b/a) ~ 0.91 and (c/a) ~ 0.84. Rotation plays a negligible role in 
the support of the accreted stellar component, which is character- 
ized by a strong radial anisotropy in its velocity distribution. This 
anisotropy grows from the inside out, from cr^ ~ 2 at in the lumi- 
nous outskirts of the galaxy to cr^ 5 near the virial radius. 

• The accreted stellar component is distinct from the dark mat- 
ter halo as well as from the satellite population, which are typi- 
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Figure 12. a-top: Radial velocity dispersion of Milky Way stellar halo trac- 
ers (filled circles, from Battaglia et al 2005). These data are approximated 
with a simple model where ct^ oc (1 + r/ro)"^, shown by the thin line 
through the data points. The three other curves (scale at right) illustrate the 
radial dependence of the right-hand side terms in Jeans' equation (eq.|2j, de- 
rived assuming a Sersic profile for the tracers, h-bottom: Solid squares show 
the total circular velocity of the Milky Way derived from the top panel. The 
top solid curve shows a fit assuming an (adiabatically contracted) NFW 
dark matter halo and an exponential disk. The contribution of the dark mat- 
ter halo (before and after contraction), as well as that of the disk, are also 
shown. See text for further details. 

cally less concentrated and more isotropic. Most stars in the outer 
halo formed in progenitors that have since merged with the central 
galaxy. Only a small fraction of the outer halo stars are contributed 
through "harassment" of satellites surviving as self-bound entities 
until the present. 

These properties are in broad agreement with recent observa- 
tions of the outskirts of spiral galaxies as well as of giant stars in 
the outer regions of the Milky Way and of M3 1 . In particular, they 
show that the outer light "excess" over extrapolations of the inner 
luminous profile recently reported for M31 and other spirals is a 
generic feature of galaxies formed hierarchically, and that the de- 
cline in velocity dispersion seen in the outer Milky Way halo tracers 
is consistent with the presence of an extended, massive halo of dark 
matter (albeit one of perhaps lower mass than commonly assumed). 
They also imply that the intergalactic stellar component of galaxy 
clusters may be more intimately related to the central galaxy than 
to the galaxy cluster population as a whole. These results illustrate 
ways to unravel the clues to the tumultuous accretion history of 
individual galaxies contained in the stars ejected beyond their bor- 
ders. 
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Table 1. Structural parameters of dark and luminous components of simulated galaxies 



Label 


^vir 

kpc 


[10" Mq] 


[10" Mq] 


iWstr 

[10" Mq] 


[10" Mq] 


[10" Mq] 


-Mout 
[10" Mq] 




j^acc 
J out 


KIAl 


391.53 


29.47 


25.96 


3.32 


2.36 


0.47 


0.48 


0.38 


0.91 


KIA2 


266.41 


9.28 


7.88 


1.29 


1.04 


0.15 


0.10 


0.30 


0.95 


KIA3 


267.57 


9.41 


7.92 


1.25 


1.02 


0.10 


0.12 


0.42 


0.96 


KIA4 


350.05 


21.06 


17.96 


2.75 


1.91 


0.51 


0.33 


0.56 


0.94 


KIA5 


316.30 


17.61 


15.13 


2.25 


1.68 


0.35 


0.23 


0.40 


0.96 


KlBl 


394.00 


31.53 


27.19 


3.53 


2.60 


0.55 


0.37 


0.39 


0.93 


KIB2 


269.59 


9.62 


8.18 


1.33 


1.04 


0.17 


0.11 


0.38 


0.97 


KIB3 


267.83 


9.43 


7.94 


1.23 


0.97 


0.13 


0.12 


0.49 


0.95 



Mdi-k and Mstr are, respectively, the dark and stellar mass within the virial radius, rvjr, defined to encompass a region 100 times denser than the critical 
density for closure. is the total mass of stars within the luminous radius, rim,j = 20 kpc. Msat is the mass of stars in satellites. Mout is the mass of 

stars outside the luminous radius but not in satellites, /gj"^ and /^^Jr are the fractions of accreted stars in the inner (r < r\^^) and outer (r > r\^^) galaxy, 

respectively. 

Table 2. Parameters of Sersic-law fits to i?-band surface brightness profile of accreted stars in the outer (30 < i?/kpc < 130) galaxy. 



Label 


n 


kpc 


Meff 

mag/asec^ 


rfit 
-'"tot 

[lOiOL©] 


r acc 
-'"tot 

[IOIOLq] 


KIAl 


5.01 


UAl 


23.62 


3.15 


3.68 


K1A2 


6.13 


6.05 


23.48 


0.85 


1.03 


K1A3 


5.73 


5.45 


22.93 


1.10 


1.44 


KIA4 


5.99 


6.33 


22.36 


2.61 


3.52 


KIA5 


8.27 


11.06 


24.54 


1.23 


2.07 


KIBl 


6.95 


8.49 


23.10 


2.57 


3.78 


KIB2 


6.10 


6.07 


23.30 


1.01 


1.32 


KIB3 


5.96 


5.62 


23.15 


0.97 


1.90 



is the total luminosity of the Sersic fit to the outer profile. L^°^ is the total luminosity of the accreted stellar component. 
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